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LONG-TERM  GOALS 

The  goal  of  this  research  is  to  develop  greater  understanding  of  the  dynamics  of  fluid  mud  and  its  role 
in  the  transport  and  deposition  of  sediment  in  coastal  environments.  In  particular,  we  seek  greater 
understanding  of  the  processes  that  influence  the  formation  and  maintenance  of  fluid  mud  in  energetic 
environments. 

OBJECTIVES 

The  research  is  a  process-based  study  that  addresses  three  primary  objectives: 

•  Determine  controlling  factors  in  the  formation  and  destruction  of  fluid  mud  under  a  sheared 
flow 

•  Verify  Richardson  number  dependence  for  suppression  of  turbulence  and  carrying  capacity 
of  a  high-concentration  suspension 

•  Evaluate  effects  of  mixed  grain  size  on  high-concentration  suspensions 


APPROACH 

To  evaluate  the  controls  on  fluid-mud  fonnation  and  the  influence  of  sediment-induced  stratification 
on  flow,  a  suite  of  instrumentation  has  been  assembled  to  make  in  situ  measurements  of  vertical 
gradients  of  velocity,  suspended-sediment  concentration  (SSC),  and  fluid  density  throughout  the  water 
column.  In  addition  to  these  parameters,  the  thickness  of  a  fluid  mud  layer  through  accelerating  and 
decelerating  flows  has  been  measured  in  an  environment  characterized  by  high  concentration  bottom 
suspensions. 
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The  Petitcodiac  River  located  in  the  Upper  Bay  of  Fundy  was  selected  as  the  study  site.  The 
Petitcodiac  Estuary  is  a  macro-tidal  environment  that  has  been  modified  by  the  construction  of  a 
causeway  in  the  late  1960s.  During  construction,  depositional  rates  downstream  of  the  causeway  were 
on  the  order  of  1  cm  day'1  resulting  in  a  decrease  in  the  cross-sectional  area  of  up  to  90%  within  one 
year  of  closure.  Typical  maximum  tidal  currents  are  on  the  order  1-2  m  s'1  and  suspended-sediment 
concentrations  regularly  exceed  10  g  T1  and  can  reach  300  g  l'1  (Curran  et  ah,  2004).  At  slack  tide, 
sediment  settles  rapidly  forming  fluid  mud  layers  on  the  order  of  1-2  m  thick.  In  April  2010,  the 
causeway  gates  were  permanently  opened  to  re-establish  tidal  flow  throughout  the  estuary.  The  effect 
of  the  increase  in  tidal  flow  on  the  formation  and  resuspension  of  the  dense  bottom  layer  will  be 
examined  in  November  2010. 

The  Petitcodiac  River  near  Moncton,  New  Brunswick  serves  as  an  ideal  natural  lab  where  the 
formation  and  destruction  of  fluid  mud  can  be  studied  under  a  sheared  flow  with  a  quasi-steady  current 
and  without  the  complicating  effect  of  waves.  Measurements  are  carried  out  from  a  bridge  that 
provides  a  solid  framework  for  instrument  deployment  thus  eliminating  many  of  the  challenges  of 
shipboard  operations  in  high-current  environments.  The  Petitcodiac  has  a  distinct  advantage  over 
laboratory  studies  in  terms  of  scale  and  avoids  the  practical  issues  of  dealing  with  the  quantities  of  mud 
required  to  create  a  fluid  mud  layer  and  maintaining  flume  equipment  at  these  high  concentrations  of 
suspended  sediment. 


Figure  1:  Photographs  of  the  modified  SUBS  packages  deployed  in  the  Petitcodiac  River.  Left  panel  (A)  shows  the 
interior  of  SUBS  I  with  the  RBR  CTD  and  sample  bags  exposed.  The  Marsh  McBirney  electromagnetic  current  meters 
are  located  on  the  front  of  the  package  near  the  railing  and  the  OBS  and  sample  port  are  situated  on  the  opposite  side 
near  the  front.  The  right  panel  (B)  shows  a  SUBS  modified  to  hold  a  Sontek  AD  V  at  the  front,  a  D&A  OBS5  in  the 
center  and  a  downward  looking  Nortek  AquaDopp  at  the  tail  (SUBS  II).  Also  seen  in  this  panel  is  a  surfboard  with  a  1.2 
MHz  RDI ADCP  and  dual  frequency  Knudsen  sounder. 

For  water-column  sampling  in  high  flow  conditions,  a  streamlined  underwater  buoyancy  (SUBS) 
package  designed  at  the  Bedford  Institute  of  Oceanography  was  modified  for  use  as  an  instrument 
profiler  (Fig.  1  A;  SUBS  I).  The  addition  of  lead  weights  to  the  package  in  place  of  the  nonnal 
buoyancy  permits  vertical  profiling  from  a  bridge  at  current  speeds  in  excess  of  1  m  s'1.  A  CTD  and  an 
optical  backscatterance  sensor  (OBS)  are  mounted  through  the  front  and  side  of  the  package 
respectively.  A  sample  port  is  co-located  with  the  OBS  and  connected  to  a  series  of  pumps  to  provide 
suspended  sediment  for  OBS  calibration  and  grain  size  analysis.  Current  shear  was  initially  determined 
using  two  Marsh  McBirney  electromagnetic  current  meters  (EMCMs)  mounted  60  cm  apart  on  the 
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front  of  the  CTD  sampling  package.  In  2009  these  measurements  were  augmented  with  the  addition  of 
a  second  SUBS  package  that  was  modified  to  accept  a  Sontek  ADV,  a  D&A  OBS5,  and  a  Nortek 
Aquadopp  profiler  (Fig.  IB;  SUBS  II).  To  observe  the  formation  of  the  lutocline  and  continuously 
monitor  flow  velocities,  a  dual  frequency  Knudsen  echo  sounder  and  an  1 .2  MHz  RD  Instruments 
ADCP  were  mounted  on  a  surfboard  tethered  to  the  bridge  next  to  the  location  of  the  profiling 
package.  In  2009  a  second  1.2  MHz  ADCP  and  a  2.5  MHz  rotary  sonar  were  deployed  using  a  second 
surfboard  (Fig.  2).  The  ADCPs  and  Knudsen  echo  sounder  sample  continuously,  SUBS  I 
(CTD/OBS/EMCMs/in  situ  samples)  was  deployed  approximately  every  30  minutes,  and  SUBS  II 
(ADV/Aquadopp/OBS5)  was  deployed  nearly  continuously  (every  15  minutes). 


Figure  2:  Photo  of  the  instrumented  surfboards  deployed  from  the  bridge  in  the  Petitcodiac  River.  Upper  board  consists 
of  a  1.2  MHz  RDI  ADCP  and  a  Dual  Frequency  (50  and  200  kHz)  Knudsen  sounder.  The  lower  board  has  a  2.5  MHz 
rotary  sonar  and  1.2  MHz  RDI  ADCP  mounted  on  it.  The  wake  between  the  boards  is  caused  by  a  SUBS  package  just 
below  the  surface. 

All  work  is  being  conducted  collaboratively  between  Tim  Milligan  and  Brent  Law  of  the  Bedford 
Institute  of  Oceanography  (BIO),  Gail  Kineke  of  Boston  College  (BC)  and  Alex  Hay  (Dalhousie). 
Technical  support  in  2009  was  provided  by  BIO  Co-op  students  Vanessa  Page  and  Karen  Devitt, 
Boston  College  students  Michele  Lennon  (Masters)  and  Christie  Hegermiller  (undergraduate),  and 
Dalhousie  graduate  student  Richard  Cheek 

WORK  COMPLETED 

The  original  SUBS  profiling  package  was  redesigned  in  2008  to  accept  an  RBR  XRX-620  CTD  and 
data  logger,  and  a  new  control  system  was  developed  for  the  pump  samplers.  A  second  SUBS  was 
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modified  to  accept  a  Sontek  ADV,  a  D&A  OBS5  through  the  bottom  of  the  package  and  a  downward 
looking  Nortek  AquaDopp  profiler  was  mounted  on  the  tail. 

An  experiment  was  conducted  from  June  17-22,  2009.  From  June  17-21,  conditions  were  dominated 
by  tidal  flow  with  very  low  fresh  water  input.  On  June  22  fresh  water  was  released  from  the  control 
structure  as  part  of  a  fish  migration  strategy.  Conditions  after  June  22  were  dominated  by  fresh  water 
flow  with  salt  water  intrusion  only  occurring  towards  the  end  of  the  flood. 

Initial  casts  with  the  CTD/OBS  SUBS  package  were  made  just  prior  to  the  passage  of  the  tidal  bore 
and  continued  approximately  every  30  minutes  until  water  levels  during  the  ebb  were  less  than  1  m. 

The  ADV/OBS5/AquaDopp  SUBS  was  profiled  nearly  continuously  over  the  tidal  cycle  with  1  minute 
integrated  samples  being  collected  approximately  every  0.5  m. 

Continuous  observations  of  current  speed  and  depth  of  the  lutocline  were  made  using  the  surfboard 
mounted  ADCPs  and  the  dual  frequency  Knudsen  sounder.  The  combined  measurements  of  the  SUBS 
I  package  and  the  Knudsen  profiler  are  summarized  in  Figure  3.  The  close  comparison  of  water  depth 
from  the  two  different  instrument  packages  ~  10  m  apart,  as  well  as  good  agreement  in  the  depth  to  the 
lutocline  using  acoustic  (dual  frequency  echo  sounder)  and  optical  (optical  backscatterance)  methods  is 
encouraging  for  combining  measurements  from  different  instruments  to  obtain  coverage  over  a  range 
of  conditions.  The  acoustic  instruments  are  helpful  at  defining  gradients  in  the  water  column,  e.g.  the 
lutocline,  but  are  limited  by  high  concentrations  for  actual  concentrations.  The  ADV  potentially  can 
provide  turbulence  measurements,  however  lab  experiments  suggest  only  to  SSCs  of  ~  50  g  l'1 
(Velasco  and  Huhta,  2008)  and  our  field  observations  suggest  to  ~30  g  1 1  or  less. 

Calibration  of  the  OBS  is  accomplished  using  the  in  situ  sediment  samples  that  range  from  170  mg  T1 
to  nearly  450  g  l'1.  This  calibration  takes  advantage  of  the  non  linear  response  of  the  OBS  at 
concentrations  greater  than  ~5  g  l1  that  was  demonstrated  for  fluid  mud  on  the  Amazon  Shelf  (Kineke 
and  Sternberg,  1995).  Disaggregated  Inorganic  Grain  Size  (DIGS)  analysis  for  the  suspended  sediment 
samples  was  completed. 

RESULTS 

Conditions  observed  during  the  June  2009  experiment  were  similar  to  those  in  2007.  Initially  the 
system  was  dominated  by  tidal  flow  which  resulted  in  salinities  ranging  from  1  to  14  psu  and  near 
bottom  SSC  values  in  excess  of  400  g  U .  The  passage  of  the  tidal  bore  created  mixed  water  column 
with  SSC  concentrations  on  the  order  of  30-100  g  l1  (Fig.  4).  As  current  velocity  decreases,  a  lutocline 
forms  and  near  bottom  concentrations  increase  rapidly  to  greater  than  100  g  l1  (Figs.  3,  4).  The  the 
initial  settling  rate  during  the  formation  of  the  lutocline  is  ~6  cm  s'1  .  The  fluid  mud  layer  remains 
coherent  throughout  the  ebb.  After  the  release  of  water  from  the  head  pond  on  June  22,  salinities  were 
below  1  psu  and  maximum  near  bottom  SSC  values  decreased  to  52  g  l'1. 

Quality  control  analysis  is  currently  underway  to  compare  shear  measurements  from  paired  EMCMs, 
ADV  and  ADCP  mean  profiles.  Some  differences  between  velocity  results  are  under  scrutiny  in  order 
to  assimilate  shear  measurements  with  prior  experiments,  and  complete  calculations  of  gradient 
Richardson  numbers.  To  date,  the  threshold  condition  suggested  by  Trowbridge  and  Kineke  (1994)  for 
the  carrying  capacity  of  a  flow  appears  to  hold,  but  implies  a  critical  gradient  Richardson  number 
between  0.25  and  1  for  the  rapid  fonnation  of  a  high  concentration  layer  (Heath  2009). 
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Petitcodiac  June  21  2009 


Figure  3:  Summary  of  sampling  from  the  SUBSI  profiler  and  Knudsen  dual  frequency  echo  sounder  mounted  on  a 
surfboard  10  m  away.  The  close  comparison  of  water  depth  from  the  two  different  instrument  packages,  as  well  as  good 
agreement  in  the  depth  to  the  lutocline  using  acoustic  (dual  frequency  echo  sounder)  and  optical  (optical 
backscatterance)  methods  is  encouraging  for  combining  measurements  from  different  instruments  to  obtain  coverage 
over  a  range  of  conditions. 


Figure  4:  Pumped  suspended  samples  were  used  to  convert  optical  backscatterance  measurements  to  suspended 
sediment  concentration,  although  the  calibration  result  was  noisy,  with  large  uncertainty  at  high  concentrations.  Mean 
velocity  profiles  obtained  by  the  AD  V  profiling  package  (SUBSII)  show  strong  flood  currents,  ~  1  m/s  at  the  surface  at 
1345,  followed  by  rapid  clearing  of  the  upper  water  column  and  formation  of  a  high  concentration  bottom  layer  by  1415. 
Currents  reverse  during  the  ebb  at  the  bottom  (beginning  at  1415)  yet  unusually  strong  winds  blowing  up  the  channel 
delayed  the  change  in  direction  at  the  surface. 
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The  evolution  of  the  DIGS  during  the  period  when  the  flood  slows  and  reverses  shows  a  dramatic 
difference  between  the  normal  flood  current  conditions  with  the  causeway  closed  June  2 1  and  the 
freshwater  release  event  of  June  22  (Fig. 5).  While  the  concentration  immediately  after  the  passage  of 
the  tidal  bore  varied  by  an  order  of  magnitude  (solid  black  lines  in  Fig.  5),  the  size  distribution  of  the 
material  suspended  by  flooding  current  was  similar.  Salinity  on  June  22  was  on  the  order  of  20  PSU 
but  was  <1  PSU  during  the  discharge  of  river  water.  Key  differences  in  settling  behaviour  following 
resuspension  and  mixing  by  the  flood  tide  are  the  lack  of  floe  settling  on  the  22nd  shown  by  the 
approximately  constant  DIGS  distribution  and  concentration  over  the  period  of  low  flow  speed  and  the 
evolution  of  the  modal  peak  in  the  distributions  for  the  flocculated  samples  collected  at  14:15,  14:30 
and  15:00  on  the  2 1st  at  10  and  50  cm  above  bottom  (blue  and  green  curves,  respectively).  In  the 
flocculated  suspension  below  the  lutocline,  the  concentration  in  the  <20  pm  size  classes  remains 
constant  whereas  the  concentration  in  the  >20pm  fraction  increases.  This  would  suggest  that  within 
the  dense  bottom  layer  floe  settling  is  completely  hindered  but  single  grains  can  penetrate  to  the 
bottom,  although  their  settling  is  also  hindered.  In  the  unflocculated  suspension,  there  is  sufficient 
time  over  the  45  minute  sampling  period  for  all  particles  >20  pm  to  be  removed  from  suspension.  The 
effective  gelling  of  the  floes  below  the  lutocline  combined  with  the  continued  settling  of  the  single 
grains  through  the  dense  bottom  layer  could  explain  the  formation  of  the  varves  observed  in  the 
Petitcodiac  Estuary  (Fig.  6).  These  varves  separate  along  a  coarse  silt  layer  which  abruptly  changes  to 
fine  grained,  unsorted  mud. 


June  21 , 2009  June  22,  2009 


Figure  5:  Disaggregated  inorganic  grain  size  distributions  for  samples  collected  at  30  cm  below  the  surface  (red),  10  cm 
above  bottom  (  blue)  and50  cm  above  bottom(green)  during  maximum  flood  when  flow  speeds  went  from  approximately 
-0.25  (solid  line)  to  0  (dashed  line)  to  +0.25  cm  s'1  (dashed  and  dotted  line).  The  black  line  is  the  size  distribution  of  the 
suspension  immediately  after  the  passage  of  the  tidal  bore.  The  surface  salinity  range  on  June  21  was  14-19  and  on 
June  22was  <1.  Data  for  the  50  cm  above  bottom  samples  on  June  22  were  omitted  for  clarity  but  were  the  same  as  those 
for  10  cm  above  bottom. 


6 


Figure  6:  Image  of  the  varved  sediment  that  accumulates  on  the  banks  of  the  Petitcodiac  River  at  Moncton  during 
periods  of  rapid  accumulation.  The  van’es  separate  along  a  coarse  silt  layer. 

IMPACT/APPLICATION 

Our  recent  observations  are  refining  our  understanding  of  current  shear  at  the  lutocline  surface.  These 
observations  improve  our  understanding  of  the  formation  and  maintenance  of  fluid  mud  in  sheared 
flows  and  to  further  test  threshold  conditions  for  the  suppression  of  turbulence  and  carrying  capacity  of 
turbulent  flows  (Trowbridge  and  Kineke,  1994).  Multiple  measurements  with  different  instruments 
provide  a  unique  dataset  on  how  acoustic  instruments  behave  in  high  concentration  sediment 
suspensions.  Our  data  confirm  the  importance  of  salinity  and  sediment  concentration  in  the  formation 
of  fluid  mud  layers  and  provide  a  basis  for  investigating  the  dynamics  of  these  layers  in  macrotidal 
estuaries. 

RELATED  PROJECTS 

G.  Kineke  is  a  co-PI  on  an  ONR  MURI  project,  “Mechanisms  of  Fluid-Mud  Interactions  Under 
Waves.”  This  project  includes  investigators  from  Johns  Hopkins  (Dalrymple,  Shen),  Woods  Hole 
Oceanographic  Institution  (Trowbridge,  Traykovski),  MIT  (Liu,  Mei,Yue)  and  Memorial  University 
(Bentley).  The  major  objective  of  this  study  is  to  examine  the  various  mechanisms  of  water  wave 
dissipation  over  muds  and  field,  laboratory,  and  theoretical  approaches  will  be  employed. 
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T.  Milligan  is  a  co-PI  on  the  ONR  Tidal  Flats  project  in  Willapa  Bay  looking  at  sediment  flux  and  its 

impact  on  sediment  strength. 

REFERENCES 

Curran,  K.J.,  T.G.  Milligan,  G.  Bugden,  B.  Law,  and  M.  Scotney.  (2004).  Suspended  Sediment,  Water 
Quality,  and  Hydrodynamics  of  the  Petitcodiac  River  Estuary,  New  Brunswick  (2002-2003). 

Can.  Tech.  Rep.  Fish.  Aquat.  Sci.  2516:  xi  +  88  p. 

Heath,  K.M.  (2009)  Fluid  Mud  Fonnation  in  the  Petitcodiac  River,  New  Brunswick,  Canada.  Masters 
Thesis,  Boston  College,  Chestnut  Hill,  MA.  125  pp. 

Kineke,  G.C.  and  R.W.  Sternberg  (1992).  Measurements  of  high  concentration  suspended  sediments 
using  the  Optical  Backscatterance  Sensor.  Marine  Geology,  108:253-258. 

Trowbridge,  J.H.  and  G.C.  Kineke  (1994).  Structure  and  dynamics  of  fluid  muds  on  the  Amazon 
continental  shelf.  Journal  of  Geophysical  Research,  99(Cl):865-874. 

Velasco,  D  and  Huhta,  C.  (2008).  ADV  performance  in  fluid  mud.  Poster  presentation,  Biennial 
Conference,  Physics  of  Estuaries  and  Coastal  Seas,  Liverpool. 


